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Abstract. We consider a nonlinear structured population model with a distributed 
recruitment term. The question of the existence of non-trivial steady states can be 
treated (at least!) in three different ways. One approach is to study spectral prop- 
erties of a parametrized family of unbounded operators. The alternative approach, 
on which we focus here, is based on the reformulation of the problem as an integral 
equation. In this context we introduce a density dependent net reproduction rate 
and discuss its relationship to a biologically meaningful quantity. Finally, we briefly 
discuss a third approach, which is based on the finite rank approximation of the 
recruitment operator. 



1. Introduction 

In mathematical epidemiology the basic reproductive number, often denoted by Rq, is 
the expected number of secondary infections that a single infected individual will cause in 
a completely susceptible population when there is no control or intervention, see e.g. [25| . 
In the mathematical context these ideal assumptions mean basically that the underlying 
model is a linear one (or indeed it is a linear approximation of a nonlinear one) . Hence it 
is expected that the infection will persist if i?o > 1- It is also expected that larger values 
of Rq will result in a major epidemic and significant intervention efforts will be needed to 
control the spread of the disease. In a mathematical model Rq is in fact often introduced 
as a threshold parameter that allows one to determine the stability of the disease free 
steady state, the spectral radius of a bounded positive linear operator, see 

In population ecology/dynamics, the value Rq is often related to the measure of re- 
productive success of an individual in an ideal environment, see e.g. [llj. It can be then 
used as a bifurcation parameter to study the existence and local asymptotic stability (in 
some cases global asymptotic stability) of the extinction (or in some case the positive) 
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steady state. In the context of the fohowing hnear age-structured population model 

Pt{a,t) +pa{a,t) = -fi{a)p{a,t), p(0,i) = j (3{a)p{a,t)da, (1.1) 
one introduces the quantity 

R ~ /3(a) exp < — / fi{7') dr > da, 



which is the expected number of newborns of an individual to be produced in her lifetime. 
It is then shown that p{a, t) ~ e''*p*(a), as t — >■ oo, where r and R — 1 have the same sign 
(see (i^l for general results of this type) . This asymptotic property of solutions is called 
asynchronous exponential growth and it can be effectively characterized in the framework 
of semigroup theory, see e.g. [H, [T^iE^- We also note that equivalently R can be defined 
as the spectral radius of the Volterra integral operator 

B{t)^ y" /3(a) cxp /i(r) drj B(t - a) da, (1.2) 

defined by the right hand side of the renewal equation (|1.2p . where B{t) = p{0,t). The 
advantage of formulating the linear age-structured problem (|l.ip as a renewal equation 
is that integral operators are usually nicely behaved (for example bounded) in contrast 
to unbounded differential operators. 

As the above simple linear example illustrates, net reproduction numbers are expected 
to play a key role in the analysis of both discrete and continuous population models. 
Recently, this has become a resurgent topic of interest. In the context of structured 
epidemiological models we refer the reader to the recent papers by Bacaer at al. [3, Q , 
while in the context of discrete structured population models we refer to |llMl3| . It is 
clear that if one adds nonlinearities to the simple age-structured model (jl.ip . it cannot 
be expected that a simple constant will determine the asymptotic behaviour of solutions, 
for example whether the population will persist or die out, in general. However, it is 
still true that questions of existence and local asymptotic stability properties of steady 
states can be naturally related to appropriately defined net reproduction functions (or 
functionals) see for example [13, [3l • 

Recenty in [isj we considered a Gurtin-MacCamy type model with a distributed re- 
cruitment term 

4^5, t) + -^ (7(s, P{t))p{s, t)) = -/i(s, P{t))p{s, t)+ [ (3{s, y, P{t))p{y, t) dy, (1.3) 
at ds Jq 

7(0,P(t))p(0,t)=0, (1.4) 

p(s,0) =po(s), P{t)^ f p{s,t)ds. (1.5) 
Jo 

Here we set the maximal size to be 1 for mathematical convenience. In this model it 
is assumed that individuals may have different sizes at birth and therefore /3(s, y, ■ ) 
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denotes the rate at which individuals of size y give rise to individuals of size s. Hence 
the non-local integral term in (|1.3p represents reproduction of the population without 
external driving of the population through immigration. Existence of solutions of the 
model above (and in fact of a much more general model) was treated in [§]. In we 
established some sufficient conditions for the existence of a positive steady state of this 
model. For a fixed P S (0, cxd) we defined a linear operator Ap hy 

Apu = --^ (7(-, P)u) - fii; P)u + / /?(•, y, P)u{y) dy, 
OS Jo 

Dom(^p) ={ue W^'\0, 1) I u{0) 0} . (1.6) 

Then we established conditions which guarantee that there exists a such that Ap^ has 
eigenvalue with a corresponding unique positive eigenvector. (See also where this 
approach was employed in case of a cyclin-structured cell population model.) To this 
end, we established that is the generator of a positive and irreducible (under some 



mild condition on /3) semigroup. In jl8j we also established that the semigroup generated 
by Ap is eventually compact, which implies that the Spectral Mapping Theorem holds 
true, and the spectrum of may contain only isolated eigenvalues of finite algebraic 
multiplicity (see e.g. [3). It then follows that the spectral bound is a dominant (real) 
eigenvalue Ap = s (Ap) of geometric multiplicity one with a corresponding positive 
eigenvector [13, Chapter 9]. Finally, we established conditions which imply that there 
exist a P+ S (0, oo) such that the spectral bound s{Ap+) is negative and therefore 
the dominant eigenvalue Ap+ ~ s{Ap+) is also negative; and a P~ S (0, oo) such that 
this dominant eigenvalue Ap- = s{Ap-) is positive. Then it follows from standard 
perturbation results on eigenvalues (see e.g. (22l |) that there exists a zero eigenvalue. 

The reformulation of the steady state problem as an eigenvalue problem for a family 
of unbounded linear operators allowed us to obtain sufficient conditions for the existence 
of positive steady states. Unfortunately, the biological motivation is somewhat lost in 
the approach above. In the present paper we treat the question of existence of positive 
steady states of model (jl.3p - (|1.5p using a different method. Our motivation is to establish 
the existence of positive steady states via an appropriate net reproduction function as 
in the case of the basic Gurtin-MacCamy model with one state at birth. Nevertheless, 
since local (and in some cases global!) asymptotic stability of equilibria can naturally be 



related to an appropriate net reproduction function, see e.g. jl7l-ll 



2. Existence of steady states and the net reproduction function 

In this section we discuss the existence of steady states of model (|1.3p - (|1.5p by first 
(re) formulating it as an integral equation. To this end we introduce the function 

B*{s):=B{p,,s)^ ( f3{s,y,PMy)dy. (2.7) 
Jo 
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With this notation an imphcit solution of the steady state equation can be obtained as 
= Texp/- r t^^flfA dAB*{r)dr, (2.8) 

1{V,P*) Jo I Jr l{x,P^) J 

which yields the integral equation 

B*{s)= C^-^^f^ fcxpl- r/f(^d4s*(0drdy, .e[0,l]. (2.9) 

Jo l{y,P*) Jo { Jr 1{X,P*) J 

Using the notation 

K{s,r,P.,)= cxpi- — — ^dx)- — - — B^'^y^ 

Jr { Jr 1{X,P*) J l{y,P*) 

and a change of variables in the integration, we can recast equation (j2.9l) in the more 
economic form: 



B*{s)^ I K(s,r,P^)B*(r)dr, .s6[0,l]. 
Jo 



(2.10) 



The existence of a non-trivial solution of equation (j2.10p is necessary for the existence of 
a non-trivial steady state . On the other hand, once a non-trivial solution of p.lOp is 
found we can substitute this directly into (2.5) to get a positive steady state. Also note 
that if the integral equation (|2.10p has a non-negative solution then any positive scalar 
multiple of that is a solution, hence we can normalize B* such that the steady state will 
also satisfy the following condition 



P* = p*iy) dy- 

Jo 



We note that in the spirit of the approach we employed in [18[ , problem (|2.10p can be 
treated as an eigenvalue problem, now for a bounded operator. In particular for a fixed 
"environment" P we consider the linear integral operator 



Cpx= K{-,r,P)x{r)dr, for xeX. (2.11) 
Jo 



More precisely we consider a family of integral operators parameterized by P. The exis- 
tence of a positive (not necessarily strictly positive, only in the positive cone) steady state 
requires that for some P* > the operator Cp^ has eigenvalue one with a corresponding 
positive (not necessarily strictly positive!) eigenvector x. 

We make the following regularity assumptions on the model ingredients which are 
needed in what follows. 

7 e C2([0, 1] X [0, oo)), n e Ci([0, 1] X [0, oo)), (3 S C\[0, 1] x [0, 1] x [0, oo)), 

< 70 < 7 < r, 0< fi<M, 0< f3<B. (2.12) 

Some of the above conditions are biologically relevant, such as the assumption of non- 
negative and bounded vital rates, the others, such as differentiability, are necessary to 
discuss linear stability of steady states, see later in Section 5. 
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The relationship between problems (jl.6[) and (|2.11|) is established using the main result 
of 27[, which we recall here for the reader's convenience. 

Theorem 2.1. f f2^ . Theorem 3.5]) Let B be a resolvent-positive operator in X , s{B) < 0, 
and A = B + C a positive perturbation of B. If A is resolvent-positive then s{A) has the 
same sign as r (— CS^^) — 1. 

In our setting for every P G [0, oo) we have Ap = Bp + Cp, where is defined in 
(dH) and 

Bpu = --^{J{■,P)u)~^l{■,P)u, Dom(Bp) = {wG W^i'i(0,l)|w(0) =0), (2.13) 
OS ^ ' 

Cpu^l p{-,y,P)u{y)dy, Dom(Cp) = £^(0, 1). (2.14) 

JO 

It is clear that for every P G [0, oo) Cp is a positive operator and Bp is resolvent-positive, 
i.e. its resolvent set contains a positive halfline and its resolvent {XI— Bp)^^ is positive for 
A large enough. This is because it is shown that it generates a positive quasicontractive 
semigroup (see [l3l)- We also note that if inf{/i(s,P) | s G [0, m]} > v > Q holds then 
it is shown that Bp generates a positive contraction semigroup, in particular ujq{Bp) = 
s{Bp) < —v. We note that the assumption of a strictly positive mortality function is 
very natural from the biological point of view for most populations. Next we show that 
for every P G [0, oo) the integral operator Cp defined in (|2.1ip is in fact —CpBp^. Since 
we have 

5(r-{-r^-}"'^'- 

it is easily shown that if we define an operator as 

-1 r ( f Kx,p) 

l{;P)Jo n Jr l{x,P) 
with Dom (^Bp^) — i^(0, 1) then we have 

BpBp^u = u, (2.16) 

for every P G [0,c«), hence Bp^ is indeed the (right) inverse of Bp. If the regularity 
conditions (j2.12p hold true then it is shown that the operator Bp^ maps into the 
domain of Bp for every P. Thus, we have 

Cpu=-CvBp^u, (2.17) 

for aU u G 2.^(0, 1) and P G [0, oo). 

We also note that the problem (|2.1ip can also be treated as a "fixed-ray" problem 
in the positive cone oi X. In particular it is clear that the operator Cp (restricted to 



^p'"=rr^ / exp<j- / dx \ u{r) dr (2.15) 
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the positive cone of X) maps positive rays into positive rays, and any fixed ray yields 
an eigenvalue Xp , where Cp restricted to this fixed ray is a multiplication operator with 
Ap. The idea of discussing existence of non-trivial steady states in the framework of 
a combination of fixed point results and spectral theory seems very fruitful, sec e.g. 



Theorem 2.2. The spectrum of Cp is not empty and contains only eigenvalues of finite 
multiplicity. 

Proof. Since the kernel K is continuous for every P, the linear integral operator Cp is 
bounded for every P G [0, oo) with 



hence it is continuous and compact. 

We define a map $ : B— >-Sas<J> = 'Po£p, where B is the unit sphere intersected with 
the positive cone of X and P is a projection onto B along positive rays in X. Note that the 
set B is convex since X = is an AL-space, see e.g. We apply Schauder's Theorem 
to the map A fixed point of this map $ in turn implies the existence of a "fixed-ray" 
of the integral operator Cp. On this ray (one-dimensional subspace) the operator acts as 
a multiplication operator with a constant say Ap. This Ap is an eigenvalue of £p with 
a corresponding eigenvector which is any (non-zero) clement of the fixed-ray. □ 

Recall for example from [23| that a (non-trivial) continuous positive endomorphism O 
on the Banach lattice X is irreducible if it does not admit closed invariant ideals other 
than the trivial ones. 

Proposition 2.3. // the operator Cp is irreducible then a (Cp) = r [Cp). 

Proof. If the operator Cp is irreducible then the spectral radius is the only eigenvalue 
which has a corresponding strictly positive eigenvector. Hence any other eigenvalue in the 
spectrum of £p would have a corresponding (non-trivial) positive eigenvector vanishing 
on some set of positive measure. It is easily shown then that this element generates a 
non-trivial ideal (i.e. the minimal ideal containing the eigenvector) in the lattice X which 
is invariant under Cp. □ 
As the above proposition suggests, if the operator £p is irreducible for every P then 
it is natural to define the net reproduction rate of the standing population to be the 
spectral radius of £p, i.e. R{P) ■= r{Cp). If for example /3 is strictly positive, then 
this is clearly the case. If however /3 vanishes on some set of positive measure then Cp 
is not necessarily irreducible. In fact for example if there is a maximal offspring size M 
(less than the maximal size) then the ideal of equivalence classes of functions vanishing 
on (M, 1) will be invariant under Cp. Most importantly if there is another eigenvalue Ap 
of the operator Cp with a corresponding positive eigenvector xp (or "fixed-ray") which 
belongs to one of those ideals, then this eigenvector still yields a strictly(!) positive steady 



oo 



= max 
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state via formula (|2.8p . Hence this (or for that matter any other) eigenvalue could be 
considered as a net reproduction rate. 

The idea of proving existence of positive steady states via the net reproduction function 
works now exactly as in the case of simple Gurtin-McCamy models. Let us write: 

cr (Cp) = {riCp) ^ > Xl > ■ ■ ■ > X^p, ■ ■ ■ ; i e . (2.18) 

If Ap is any eigenvalue such that Aq > 1 and lim Ap = then it follows from the 

continuous dependence of the eigenvalue on the parameter P (see e.g. [22| ) that there 
exists a > such that Ap =1. Also Ap_ admits a positive (not strictly positive if 
i ^ 1) eigenvector which is a fixed ray for the operator Cp. We can normalize then 
the eigenvector (and write for the normalized vector) such that P» = B^{s) ds 
holds. Then formula (|2.8p yields a non-trivial steady state. 

3. Results for special types of kernels 

As we noted in the previous section, from the mathematical point of view a natural 
candidate for the net reproduction rate, at least in the case when the integral operator 
Cp is irreducible, would be the spectral radius r(Cp). From the biological point of view 
however a natural candidate is 



R{P) = [ [ K{s,r,P)drds. 
Jo Jo 



(3.19) 



This is because the kernel K gives the average density of individuals of size s produced 
by individuals who were born at size r, if the standing population is P (and would 
remain constant P). In this section we show that the definition of the net reproduction 
function discussed in the previous section (defined via eigenvalues of an integral operator) 
coincides with the biologically meaningful net reproduction rate R{P) defined in p.l9p 
at least in the case of special classes of kernels of the integral operator defined in (j2.1ip . 

First let us assume that there exist a positive constant (or more precisely a positive 
function of P) Cp such that 

Cp= I K{s,r,P)ds, Vr£[0, 1]. (3.20) 
Jo 

Theorem 3.4. // the kernel K satisfies p.20p then for every P e [0, oo) 

r{Cp)^J J K{s,r,P)drds (3.21) 

Proof. First we show that if p.20p holds then R{P) < r{Cp). To this end we are 
going to utilize the minimax principles for a bounded linear operator (and its dual) on 
a Banach space, see [2^ . Let X = L^(0, 1) with K ~ i+(0, 1) the positive cone. Then 
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X* = L°°{0, 1) and let K' = L^{0, 1) denote the positive cone of X*. Let (•, •) denote 
the natural semi-inner product defined between elements of X and X* . 

Recall from [2^ that if T is a linear endomorphism on the Banach lattice X with cone 
K then for any non-zero clement x' G X* we can define 

s^,{T) = sup {v em\{Tx,x') >v{x,x'}, xeK}. 

V 

Also recall from (23j (Lemma 3.4) that if T is a bounded positive linear endomorphism 
of X and K is closed then K' is closed, and 

s^'iT)<r{T), (3.22) 

for any x' E K' . Now for any x G K, x' £ K' and P G [0, cxd) we have 

{Cpx,x')^[ [ K{s,r,P)x{r)drx'{s)ds^ [ [ K{s,r,P)x'{s)dsx{r)dr. (3.23) 
Jo Jo Jq Jq 

In particular from equation (|3.23p for a;' = 1 G K' we have 

r-l 



{jCpx,l) > Cp / x{r)dr, 
Jo 

for any x €z K. Hence for every P G [0, oo) we have 

RiP) <si{Cp) <r{Cp). 
The following on the other hand shows that r{Cp) < R{P) holds. 

r{Cp) ^ r{C*p) < \\C*p\\^ = mux [ A'(r, s, P) dr = max Cp = P(P). (3.24) 

s6[o,il7o se[o,il 

□ 

Remark 3.5 From the biological point of view condition p.20|) implies that individuals 
born at different sizes have the same net reproductive rate. 

Theorem 3.6. (Theorem 3.4*) If the kernel K satisfies 

C*p^ [ K{s,r,P)dr, sG[0, 1], P G [0, oo) (3.25) 

then for every P G [0, cxd) 

r{Cp)= f f K{s,r,P)drds^C*p^:R{P). (3.26) 
Jo Jo 

Proof. Similarly as in Theorem l3.4l first we show that if (|3.25p holds then R{P) < r {Cp). 
To this end we are going to utilize again a minimax principle from (23j . Recall from 
[23! that a subset H' C K' is called /v-total is and only if from the relations {x,x') > 
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for all x' S ff' it follows that x € K. Note that in our setting K' is ii'-total, since K is 
closed. For any x G K define 

r^(T) = sup {w e M I {Tx - wx) G K} . 



Lemma 3.1 in (23| states that if H' C K' is K-tota\ then 

r^(T) = sup {u e R I {Tx, x') > v{x, x'), x' G H'} . 

V 

Also recall from (Lemma 3.3) that if T is a bounded positive linear endomorphism 
of X and K is closed then for any ^ x € K we have 

r^T) < r{T). (3.27) 

Again we choose a; = 1, for which for every x' £ K' we have 

{Cpl,x') >C*p I x'{s)ds. 



Hence for every P G [0, oo) we have 

R{P)<n{Cp)<r{Cp). 

To show r{Cp) < R{P) we utilise again a result from [2^. Recall that for a positive 
bounded linear endomorphism T we have 

s^'{T) = inf {r eR\T{x,x') > {Tx,x'), x £ K} . (3.28) 

r 

Again let us choose x' = 1 and it is easy to see that Cp{x, 1) > (£p, 1) for every x £ K, 
hence we have 

r{Cp) < s\Cp) < R{P). (3.29) 

□ 

Remark 3.7 From the biological point of view condition p.25p implies that individuals 
produce the same amount of offspring of different sizes during their lifetime. 

Next we return to the case when the integral operator is irreducible, since in this case 
as we have shown earlier, the spectrum only contains the spectral radius r{Cp). Recall 
from [2^ that the integral operator £p is irreducible if and only if for every S C [0, 1] 
(where S has positive Lebesgue measure) we have 

j J K{x,y,P)dydx > 0. (3.30) 
[oa]\s s 



It is relatively simple to obtain the following result. 
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Proposition 3.8. // the kernel K satisfies (j3.30p for every P G [0, oo) then 

min / A'(x, ?/, P) dx < r(£p) < max / K{x,y)<lx, (3.31) 
V Jn y Jo 

min / K{x, y, P) dy < r{C*p) < max / K{x, y) Ay. (3.32) 
^ Jo ^ Jo 

Proof. Since Cp is irreducible, let / G i^(0, 1) denote the strictly positive eigenvector 
corresponding to the spectral radius r{Cp). We have 

/ K{x,y,P)f{y)dy^r{Cp)f{x), xg[0,1]. (3.33) 

Integration of p.33p and Fubini's Theorem yields: 

f{y) I K{x,y,P)dxdy^r{Cp) f f{x)dx. (3.34) 



From p.34p we immediately obtain p.3ip noting that the kernel K is positive. 

Next we note that r{Cp) = r{C*p), and that if K satisfies p.30p then C*p is also 
irreducible. Finally, we proceed similarly as above to obtain the inequalities p.32p for 
riC*p). □ 

Remark 3.9 Note that the inequalities p.3ip - p.32p together with the assumption of 
Theorem 13.41 and also that of Theorem 13.61 would imply trivially that r{Cp) = R{P)- 
However, Proposition l3 . 81 onlv holds for an irreducible operator and neither the condition 
of Theorem 13.41 nor of Theorem 13.61 imply that Cp is irreducible. 



Remark 3.10 It is clear that we cannot expect r{Cp) ~ K{x,y, P) dy dx to 

always hold, since it generally does not hold even for (irreducible) 2x2 matrices. 



4. Finite rank approximation of the recruitment term 

In this section we briefly outline another approach to treat the steady state problem, 
which relies on the approximation of the recruitment operator with finite rank operators. 
As we noted before the general model (|1.3p - (|1.5p cannot be solved explicitly even for a 
time independent solution. Instead we consider the approximate problem with a fertility 
function /?" defined as 

n 

/3"(s,2/,P) = 5^/3,(s)^,(y)^,(P), 

r=l 

with /3r, /3r, /3 G C". We note that for this type of separable fertility function the recruit- 
ment integral operator is of finite rank (at most n). Hence for this fertility function /3" 
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we find the solution (expUcitly) of the steady state equation as 

n 

pl'(s) = ^ QTPr{P:)Fr{s, PD, (4.35) 



where 



Fr{s,P)^ exp<^-/ — — dx^— — — dy 

Jo I Jy l{x,P) J 7[y,P) 

r r f4^, \ My) ^ 
[ Jy l{x,P) J 7(s,-P) 



and 



Multiplying equation (|4.35p by I3j{s) and integrating from to to we obtain an n- 
dimcnsional system 

where Q" = (Qi'", • • • , Q"'")"'" is an n-dimensional vector for every n and Op is an rt x n 
matrix valued function with positive elements ofj{P), where 

ol{P)=MP) Ms)F,{s,P)ds, i,j = l,...,n. 
Jo 

It can be shown that a positive stationary solution to the approximate /3" problem exists 
if and only if there exists a > such that the matrix Op has eigenvalue one with 
a corresponding strictly positive eigenvector. Note that similar to the general case we 
need to assure that, now a finite-dimensional operator, has eigenvalue 1 with a positive 
eigenvector. Hence we may call Op (purely motivated from the mathematical point of 
view) a net reproduction matrix (corresponding to the approximate /3" problem). Clearly 
0"(P) = Op is a continuous function of the variable P if /i, 7 and the /3's are continuous 
functions of P. 

Next we prove that a given bounded and measurable fertility function /3 can be ob- 
tained as a limit 

n 

f3{s,y,P)= \\mY.f^ris)My)MP) (4.36) 

r=l 

with pointwise convergence in s and y uniformly on compact intervals for P, with non- 
negative functions f3r,f3r- 

Note that if /3 is a continuous function of P then the (3k constructed below will be 
continuous functions of P for every k. To this end, for every j G N we partition the 
square [0,1]^ into 4-'"^ equal subsquares denoted by i, where 1 < fc, / < 2^"-^. Set 
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/3^^^ = /3. For each j on every subsquare Sf, ; we define recursively 
4;(P)= inf , 

k,l 

where x denotes an indicator function. Finally we re-index the sequence of functions 
XkjXi s-nd ; consecutively to obtain the functions in (|4.36|) . where we have removed 
any possible zero elements from the sequence. 

The spectral radius of the matrix Op is a monotone decreasing function of P if the 
elements ofj{P) are monotone decreasing functions of P. Hence for any fixed n one can 
establish conditions on the model ingredients that guarantee the existence of a value 
P such that the corresponding matrix Op has spectral radius 1. It follows then from 
Krein-Rutman Theory that this is an eigenvalue with a corresponding strictly positive 
eigenvector. This eigenvector then yields via formula (j4.35p a strictly positive steady 
state of the approximate problem. 

Hence from the mathematical point of view it seems reasonable for the approximate 
problem to define the net reproduction rate of the standing population as the spectral 
radius of the matrix Op. Next one may ask naturally the question whether a limiting 
linear-operator-valued function exists, that is, one would like to show that 

lim O^ = Op, 

with convergence for example with respect to the topology induced by the operator norm. 
This question is outside the scope of the present paper and is left for future work. We 
note however that this limiting operator (if it exists) admits all the nice properties as does 
the integral operator Cp. Firstly, because of the uniform convergence in P, the limiting 
operator is a continuous function of P. Op is also compact since it is a limit of operators 
of finite rank. Positivity of Op follows immediately from the fact that the positive cone 
of is closed. That is, Krcin-Rutman type theory may be again applied to study spectral 
properties of Op. The relationship however, between this parametrised family of limiting 
operators Op (acting naturally on l^), and the family of integral operators Cp (acting on 
L^) defined in (|2.1ip is far from clear. We note that to discuss the relationship between 
these two objects one would naturally need to embedd the two state spaces into a Banach 
space of measures. 

5. Concluding remarks 

In this paper we investigated a nonlinear structured population model with distributed 
states at birth. We defined a biologically meaningful density dependent net reproduction 
rate, see formula p.l9|) . We also related this net reproduction rate to a mathematically 
well grounded quantity, namely the spectral radius of a linear integral operator. This is 
motivated by the question of the existence of positive steady states, see also [3 Eli- We 



NET REPRODUCTION RATE OF POPULATIONS WITH DISTRIBUTED STATES AT BIRTH 13 



have established some conditions which guarantee that these two quantities equal. We 
also presented another approach, which is based on the approximation of the fertility rate 
in the original problem. This approximation leads then to a finite dimensional problem, 
namely to an eigenvalue problem for a matrix, which is more tractable. 



It is apparent from the recent literature, see e.g. [2|-|5|, |21|, |25| that there is a recur- 
rent interest in mathematical approaches defining net reproductive numbers for epidemic 
models. This is partly due to the number of epidemic outbreaks of a variety of diseases 
in the 21st century. We point out that most of the authors have been focusing on non- 
autonomous but mostly linear models. We also note that distributed recruitment terms, 
such as the one we considered here, naturally appear in epidemic models, see e.g. 0, 
where the structuring variable represents pathogen load. 

From the ecological or population dynamic point of view intrinsic growth rates and 
net reproductive numbers play a key role in the qualitative analysis of deterministic 
models. They are often used as bifurcation parameters when proving the appearance of 
positive steady states, see e.g. [lll - [l3l |. It is also often the case, that the value of the net 
reproduction function at 0, i.e. R{0) (or Rq) determines the local asymptotic stability 
of the extinction steady state. This is also the case at least for Gurtin-MacCamy type 
nonlinear structured population models, see e.g. [31 . There is also a possibility that if 
the net reproduction rate at the zero population density, i.e. Rq is greater than one then 
the positive steady state of the model is (in fact globally) asymptotically stable. The idea 
behind this is that since nonlinearities often represent competition effects, it is natural 
to assume that the density dependent net reproduction rate is monotone decreasing 
function of the standing population size. Unfortunately, even if this is the case, there is 
still the possibility that the positive steady state looses its stability via Hopf-bifurcation, 
and periodic (stable) cycles emerge. It is however the case, see e.g. [13 [H, [2^, that 
density dependent net reproduction functions (or functionals) play a key role in the 
local stability analysis of positive steady states of physiologically structured population 
dynamic models. 

Similarly to the the single state at birth model it is also possible to relate the value 



i?(0) to the local asymptotic stability of the extinction equilibrium. We recall from [18 1 
that the linearised system around the trivial steady state reads: 

utis,t) = -7(s,0)us(s,t) - (7s(s,0) -|-/i(s,0)) u{s,t) + / l3{s,y,0)u{y,t) dy, 



7(0,0)it(0,<) = 0, U{t)= I u{s,t)ds. (5.37) 

Jo 

The governing linear semigroup has nice spectral properties, see [l^, hence it is enough to 
consider the eigenvalue equation when determining stability conditions for the extinction 
steady state. This reads: 

-7(s,0)z;'(s)-(7s(s,0) + M(s,0)) i;(s)-t- / (i{s,y,Q)v{y)dy=^\v{s), «(0) = 0. (5.38) 
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The solution of equation (|5.38|) is found (if it exists) as 

f^ r fjy) / r ^+M(^,o)+7'(^o) ^ 1 , 

where we introduced the notation 

f{x) = / f3{x,y,0)v{y)dy. 
Jo 

Hence A G C is an eigenvalue if and only if equation (j5.39p has a nontrivial solution v. 
Mulitplying equation (|5.39p by /? and integrating from to 1 yields a similar integral 
equation to (j2.9p , the characteristic equation of the linearisation at the extinction steady 
state, which reads: 

In particular for A = equation (j5.40p is identical to problem ()2.9p whith = 0. We 
may now define a family of integral operators parametrised by A (at least for real A) A/a 
where J\f\f is defined via the right hand side of equation (j5.40p . For every A this operator 
is positive, moreover its spectral radius r{J\fx) (which is an eigenvalue) is a monotone 
decreasing function of A. We therefore conclude that at least when the kernel K (or the 
net reproduction function R) is a monotone decreasing function of P then the existence 
of a positive steady state implies that equation (|5.40p has a solution with A > hence the 
trivial steady state is unstable! Hence we showed that at least in special cases that the 
definition of R agrees with the intuitive understanding that the net reproductive number 
i?o detemines local asymptotic stability of the extinction steady state and stability may 
be lost with the emergence of positive steady state. However, we have not been able 
to prove this in the general case, thus, our future efforts will focus on the development 
of an accurate numerical method to solve this model and establish convergence of these 
numerical solutions. The numerical method will then be used to verify for more cases of 
the kernel whether such a net reproductive value being less than one results in population 
extinction and if it is larger than one than the population persists. 
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